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Abstract 

We investigate the magnetic properties of a Heisenberg ferrimagnetic multilayer by using Monte Carlo simulations. The aim of 
this work is to study the local structural anisotropy model which is a possible origin of the perpendicular magnetic anisotropy in 
transition metal/rare earth amorphous multilayers. We have considered a face centered cubic lattice where each site is occupied 
by a classical Heisenberg spin. We have introduced in our model of amorphous multilayers a small fraction of crystallized Fe-Dy 
nanoclusters with a mean anisotropy axis along the deposition direction. We show that a competition in the energy terms takes 
place between the mean uniaxial anisotropy of the Dy atoms in the nanoclusters and the random anisotropy of the Dy atoms in 
the matrix. 
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1. Introduction 

Amorphous ferrimagnetic Transition Metal/Rare Earth 
(TM/RE) multilayers such as Fe/Dy are very attractive 
since they may exhibit in some conditions interesting mag- 
netic properties: a strong uniaxial perpendicular magnetic 
anisotropy, a high coercive field and a high Curie tempera- 
ture for example. The origin of the perpendicular orienta- 
tion of the magnetisation (relatively to the plane of the sub- 
strate) is not yet clearly understood, and different models 
have been proposed to explain it. These models are based 
on an anisotropic distribution of TM-RE pairs along the 
perpendicular direction [1], dipolar interactions [2], local 
structural anisotropy [3,4] or single-ion anisotropy for com- 
positionally moduled films [5] . 

In this paper, we present a Monte Carlo (MC) investi- 
gation of a ferrimagnetic multilayer consisting in classical 
Heisenberg spins with different modulus. In the framework 
of the local structural anisotropy model, we study the influ- 
ence of the single-ion anisotropy constant on the RE atoms 
and of the crystallized nanocluster concentration on the 
magnetisation orientation of amorphous Fe/Dy multilayers 
at low temperature. 



* Corresponding author. Tel: + 33-2-32-95-51-32 ; fax: + 33-2-32- 
95-50-32 

Email address: etienne.talbotauniv-rouen.fr (E. Talbot). 



2. Model and simulation technique 

2.1. Model 

We consider a face centered cubic multilayered system 
made up of Fe and Dy atoms. The Hamiltonian Jf? of the 
system is 

<ij> i6Dy 

where Si is a classical Heisenberg spin, Jy is the nearest- 
neighbour exchange interaction, Di is the single-ion mag- 
netic anisotropy constant which occurs only for the Dy 
atoms and n.; is a unit vector along the local easy axis. The 
following magnetic parameters are those of the free atoms: 
g Fe = 2, <7Dy = 4/3 and «Sx)y = 5/2. In order to relate our 
results to real systems, the J Fe -Dy interaction is consid- 
ered to be negative, whereas </ Fc _ Fc and JDy-Dy arc posi- 
tive. As in previous studies [6] , we consider that S Fc , Jfc-Fc 
and JFe-Dy are concentration dependent while the minor 
Jny-Dy interaction is kept constant : 

S Fe (X Fe ) = 1.1- 1-125(1 -X Fe ), 

J Fc _ Fc (X Fc )/fc B = 14.86 + 86.96(1 - X Fc ) (in K), 

J Fc _ Dy (X Fc )/£:B = 1.63 - 38.41(1 - X Fc ) (in K), 

^Dy-Dy/^B = 1-16 (in K), 
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Fig. 1. Diffusion profile extracted from the experimental investigation 
of (Fe 3nm/Dy 2nm) multilayers (see ref. [8]). 

where Xp c is the Fe concentration in the Fe-Dy alloy. The 
exchange interactions have been adjusted to obtain pure 
amorphous Fc and Dy Curie temperatures. In order to de- 
scribe the amorphous nature of the multilayers, exchange 
interactions are modulated by a gaussian distribution as a 
consequence of the distribution of the interatomic distance 
between nearest neighbours [7] . Since the experimental re- 
ported values for the single-ion anisotropy constant of Dy 
may vary with two order of magnitude, this parameter will 
be considered as free in our simulations. We apply periodic 
boundary conditions in the plane and free boundary condi- 
tions in the deposition direction (arbitrarily chosen as the 
z axis) to take into account the free surfaces in real sys- 
tems. The diffusion profile along the perpendicular axis has 
been obtained from tomographic atom probe analysis on 
(Fc 3nm/Dy 2nm) multilayers [8] elaborated at 570K which 
give rise to the larger perpendicular magnetic anisotropy. 
This profile is characterized by high concentrated Fc layers 
and Dy 60 — Fe4o alloy layers with diffuse interfaces between 
them (Fig. 1). 

2.2. Nanocluster model 

In order to explain the origin of the perpendicular mag- 
netic anisotropy in Fc/Dy multilayers, we were inspired 
of the local structural anisotropy model which describes 
the existence of small crystallized zones with a structure 
RE - TM 3 or RE - TM 5 [4]. This hypothesis has been re- 
tained since the experimental investigation of the Fe/Dy 
multilayers [8,9] has allowed to reject the other models. We 
have thus introduced in our model Fe-Dy crystallized nan- 
oclusters made up of 13 atoms (an Fe central atom and its 
12 nearest neighbours). We randomly dispatch between 2 
and 4 Dy atoms on the neighbouring sites belonging to the 
upper and lower planes of the central atom (Fig. 2). Then, 
the magnetic anisotropy axes of these Dy atoms are along 
the Fe-Dy bond directions whereas the other Dy atoms dis- 
play a random anisotropy. We have to mention here that 
only a few Dy atoms (those of the crystallized nanoclus- 
ters) are characterized by anisotropy axes distributed on 
average along the deposition direction. 

2. 3. Numerical simulations 

Numerical simulations were performed by using the 
importance sampling MC procedure at each temperature 
based on the standard Metropolis algorithm [10] based 
on random trial step [11]. Our results were obtained by a 



Fig. 2. Schematic representation of the nanoclusters with an Fe 
central atom (•) and 2 to 4 Dy atoms (o) as nearest neighbours. 
The anisotropy axes of the Dy atoms are represented by the dashed 
arrows. 
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Fig. 3. Temperature dependence of the reduced magnetisation (global 
and for each sublattice) of the multilayer without anisotropy. 

slow decrease of the temperature starting from the param- 
agnetic state. At each temperature, 5 x 10 3 Monte Carlo 
steps (MCS) have been discarded for equilibration before 
averaging the physical quantities over the following 10 
MCS. This procedure allows to reach equilibrium at each 
temperature when the Hamiltonian of the system does not 
display too much frustration, which is the case here. As 
a lot of parameters are generated through random pro- 
cedures (atomic positions and exchange interaction mod- 
ulation) the thermodynamic quantities have been finally 
averaged over 20 different configurations in order to get 
reliable numerical results. 

3. Results and discussion 

As a first step, we have plotted the thermal variation of 
the total and sublattice magnetisations of the multilayer 
(consisting in 18 planes of 800 atoms) without any mag- 
netic anisotropy term (Fig. 3). The multilayer exhibits a 
single magnetic ordering at about Tq = 330K for the whole 
sample contrary to the case of multilayers with sharp inter- 
faces between Fe and Dy layers. At low temperature, the 
total magnetisation of the multilayer is dominated by the 
Dy sublattice because of its large magnetic moment, and fi- 
nally the total reduced magnetisation value at OK is around 
0.52 which corresponds to a collinear ferrimagnetic configu- 
ration. In this simulation without magnetic anisotropy, the 
orientation of the magnetisation does not have any mean- 
ing since the Hamiltonian is invariant under any rotation. 
Then, in the next step, we have studied the influence of 
the single-ion anisotropy constant Dd y on the total mag- 
netisation as a function of the temperature (Fig. 4) for the 
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Fig. 4. Temperature dependence of the total reduced magnetisation 
with a fraction of 7% Dy atoms embedded in the nanoclusters for 
different values of the single-ion anisotropy constant Du y ■ 
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Fig. 5. Schematic representation of the low temperature magnetic 
structure when -DDy is small (a) or large (b). 

multilayer with 7% of Dy atoms embedded in the nanoclus- 
ters. The total magnetisation is divided by the ferrimag- 
netic magnetisation since it corresponds to the ground state 
without anisotropy. We observe that the small values of the 
single-ion anisotropy constant -DDy lead to a strongly fer- 
rimagnctic order at very low temperature with a magneti- 
sation almost perpendicular to the layers as it will be seen 
later. The larger values of £>Dy {D-Qy/k-Q > 20K) induce a 
visible decrease of the total magnetisation at low temper- 
ature which is the signature of the major influence of the 
random anisotropy on the Dy atoms not included in the 
clusters and which gives rise to sperimagnetic structures 
(Fig. 5). Finally, the magnetic configuration is the result of 
the competition between random magnetic anisotropy on 
Dy atoms in the matrix, uniaxial anisotropy on average on 
Dy atoms in the clusters and the exchange interactions. 

In Fig. 6, we have plotted the angle 9 between the direc- 
tion of the total magnetisation and the z-axis in function 
of the anisotropy constant for three fractions of Dy atoms 
included in the crystallized nanoclusters (3.5%, 7% and 
14%). For a given concentration, wc observe a decrease in 9 
when increases from IK to around 20K which means 
an increase in the perpendicular component of the magneti- 
sation. For larger -DDy, we can see a slightly increase in 9 
indicating a smaller contribution of M z due to the effect of 
random magnetic anisotropy of the Dy matrix atoms. We 
can also notice that the angle 9 is strongly dependent on the 
nanocluster concentration. A rather small concentration of 
nanoclusters is sufficient to induce a strong perpendicular 
magnetisation. 
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Fig. 6. Angle between the direction of the magnetisation and the 
perpendicular direction of the multilayer in function of the magnetic 
anisotropy constant Drj y for three values of the fraction of Dy atoms 
embedded in the crystallized nanoclusters at T = IK. 

4. Conclusion 

This work has allowed to investigate the competition 
between uniaxial and random anisotropy terms in Fe/Dy 
multilayers as a function of the single-ion anisotropy con- 
stant value and of the crystallized nanocluster concentra- 
tion. The Fe/Dy amorphous multilayer exhibits a major 
perpendicular magnetisation for only 7% of Dy atoms in the 
nanoclusters, and for D^y/k-B > 10K. Our results evidence 
that the crystallized Fe-Dy nanocluster concentration has 
more influence on the perpendicular magnetic anisotropy 
than the anisotropy constant value. The simulation of hys- 
teresis loops is actually in progress in order to complete the 
study of the nanocluster effect on the perpendicular mag- 
netic anisotropy and to get more quantitative results to be 
compared with experimental data. 
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